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Abstract 

Small variations of the entanglement entropy 5S and the expectation 
value of the modular Hamiltonian SE are computed holographically for cir¬ 
cular entangling curves in the boundary of AdS 4 , using gravitational pertur¬ 
bations with general boundary conditions in spherical coordinates. Agree¬ 
ment with the first law of thermodynamics, SS = 5E, requires that the line 
element of the entangling curve remains constant. In this context, we also 
find a manifestation of electric-magnetic duality for the entanglement entropy 
and the corresponding modular Hamiltonian, following from the holographic 
energy-momentum/Cotton tensor duality. 
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1. Introduction 

It has been known for a long time that there is a striking similarity be¬ 
tween black hole physics and thermodynamics M, which suggests a deep 
connection between gravity and thermodynamics. This led to several at¬ 
tempts to understand Einstein’s equations as effective equations emerging 
from the thermodynamics of underlying degrees of freedom, such as [4]. On 
the other hand, AdS/CFT correspondence j5h ; 7] provides a broad framework 
allowing the description of gravitational theories with AdS asymptotics in 
d + 1 dimensions as emergent from strongly coupled conformal field theories 
in d dimensions. A fair question in the AdS/CFT framework is whether the 
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similarities between gravity and thermodynamics can be explained by con¬ 
sidering Einstein’s equations as thermodynamic relations for the conformal 
held theory degrees of freedom [8j. 

More recently, it has also been suggested that the connection between 
gravity and thermodynamics should not be attributed to thermal statistics, 
but rather to quantum statistics related to quantum entanglement physics 
[SHE]. More specifically, it has been conjectured that the entanglement en¬ 
tropy, which is a measure of entanglement between subsystems of a composite 
quantum system and it is defined for a given entangling surface that sepa¬ 
rates the degrees of freedom of the corresponding conformal held theory into 
two subsystems, is directly connected to the area of an open extremal hyper¬ 
surface in the emergent asymptotically AdS geometry whose boundary is the 
entangling surface. This conjecture, which is named after Ryu-Takayanagi 
mm, provides a quantitative tool to understand how gravitational dynam¬ 
ics emerges from thermodynamics related to entanglement in the boundary 
conformal held theory. 

So far, this programme has been advanced by comparing the variation 
of entanglement entropy to the variation of the expectation value of the so 
called modular Hamiltonian for any given entangling surface mm- The 
latter can be expressed in terms of the holographic energy-momentum tensor 
when spherical entangling surfaces are taken in Poincare coordinates [IS], 
while the former is provided by the Ryu-Takayanagi formula mm- Enforc¬ 
ing the hrst law of thermodynamics for entanglement through holography, 
imposes constraints for the metric perturbations around AdS space, which to 
linear order turn out to be Einstein’s equations satisfying Dirichlet boundary 
conditions mm- In this context, it is also known that all solutions of the 
linearized Einstein’s equations with Dirichlet boundary conditions satisfy the 
hrst law of thermodynamics. 

In this paper, we work out the holographic realization of the hrst law of 
thermodynamics for small perturbations of spherical AdS^ space-time hav¬ 
ing axial symmetry and satisfying general boundary conditions. The gen¬ 
eral framework for held equations satisfying general boundary conditions in 
AdS 4 is provided in reference 125 • The entangling surfaces are now curves, 
since the boundary conformal held theory is 2 + 1 dimensional, and they 
are taken to be circular, bounding a polar cap region, so that they respect 
the axial symmetry of the bulk geometry. Our general result is that the 
hrst law of thermodynamics for entanglement is realized holographically in 
all cases, hereby extending previous works beyond Dirichlet boundary con- 
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ditions, provided that the line element of the entangling curve is inert to the 
perturbations. 

In this context, we also examine the role of gravitational electric-magnetic 
duality to entanglement physics. It is well known that small perturbations 
of maximally symmetric spaces exhibit a rank-2 generalization of electric- 
magnetic duality, interchanging the linearized Einstein equations with the 
Bianchi identities in four space-time dimensions. This symmetry was orig¬ 
inally discussed for gravitons in Minkowski space 12 a, but it was subse¬ 
quently generalized in the presence of cosmological constant by consider¬ 
ing metric perturbations of dS 4 [23] and AdS 4 space-time [21H2T| . There, 
it was also found that electric-magnetic duality in AdS 4 has a holographic 
manifestation as energy-momentum/Cotton tensor duality. These consider¬ 
ations provide the gravitational analogue of the holographic interpretation 
of electric-magnetic duality in theories with 1/(1) gauge symmetry [28], but 
their physics in the space of boundary three-dimensional conformal theories 
still remain largely unexplored. 

Duality acts on metric perturbations by interchanging their boundary 
conditions. Hence, our interest in the holographic description of gravitational 
perturbations satisfying general boundary conditions in the spirit of reference 
[29] . Extending the applications of gravitational duality to entanglement 
entropy and related issues may shed new light into this interesting subject. 

The material of this paper is organized as follows: In section 2, we present 
a review of the notions of entanglement entropy and modular Hamiltonian, 
together with their holographic description, and include various formulae 
that will be used in the computations. I 11 section 3, we discuss the general 
theory of gravitational perturbations of AdS 4 space-time and formulate the 
linearized Einstein equations as an effective Schrodinger problem, splitting 
the perturbations into two distinct classes with opposite parity. The pre¬ 
sentation is made general, encompassing arbitrary boundary conditions. In 
section 4, we compute holographically the variations of the entanglement 
entropy and the modular Hamiltonian and compare the two expressions for 
general boundary conditions. The first law of thermodynamics for entangle¬ 
ment is verified in all cases, while describing the subtleties that go into the 
computation. In section 5, we address the role of electric-magnetic duality in 
holography and study its implications for the Erst law of thermodynamics for 
entanglement. Finally, section 6 , contains our conclusions and a short dis¬ 
cussion of open problems. There are also three appendices containing various 
technical details and formulae that are used in the main text. 
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2. Entanglement Entropy and Holography 

We present a brief account of the notions of entanglement entropy and 
modular Hamiltonian together with their holographic description in terms of 
bulk space geometry. We also derive some general formulae that will be used 
later for gravitational perturbations of AdS 4 space-time satisfying general 
boundary conditions. 

2.1. First Law of Thermodynamics for Entanglement 

Consider a composite quantum system comprising of several subsystems. 
Even if the composite system lies in a pure state, with density matrix p, 
this may not be true for its subsystems, which are hereby described by a 
density matrix equal to the partial trace of p over the degrees of freedom of 
the complementary subsystem, 

p A = Tr A cp. (2.1) 

When the complementary subsystems A and A c are not entangled, the re¬ 
duced density matrix p A also describes a pure state. Entanglement between 
systems A and A c is encoded to the spectrum of the reduced density matrix 
p A under the implicit assumption that the composite system lies in a pure 
state. The entanglement entropy is defined as the von Neumann entropy 
associated to the reduced density matrix p A , 

S A ■■= -Tr (p A In p A ) . (2.2) 

The density matrix p A is Hermitian and positive semi-definite, leading to 
the definition of the corresponding modular Hamiltonian as 

Pa '■= e~ HA . (2.3) 

Then, the entanglement entropy can be rewritten in terms of the modular 
Hamiltonian as 


= -Tr {pa hi p A ) = Tr (p A H A ) = (H A ) . (2.4) 

Small variations in the pure state of the overall system or the region A 
generate variations of the density matrix, 5p A , and, thus, the entanglement 
entropy will also be perturbed. We have 

SS A = -Tr [5 (p A lnp A )] 

= -Tr (In p A 8p A ) - Tr (Sp A ) (2.5) 

= Tr (H a 6pa) = S(H a ), 
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since the trace of the density matrix is normalized to one and thus, Tr (8pa) = 
0. Thus, the variations of the entanglement entropy and the expectation value 
of the corresponding modular Hamiltonian are equal 

5S a = 8 (H a ) = 8E . ( 2 . 6 ) 

This equation is the direct analog of the first law of thermodynamics for 
entanglement physics US HZI that leads our work. 

2.2. Ryu-Takayanagi Conjecture 

The Ryu-Takayanagi conjecture [9, [ID] connects the entanglement en¬ 
tropy of a region A defined by the entangling boundary surface dA in the 
boundary field theory to the area of an extremal co-dimension two open sur¬ 
face in the bulk gravitational dual theory with boundary dA. Specifically, 
the entanglement entropy is given by 

Sa = jj- Area (.4") , (2.7) 

4(jr7v 

where H extr is the corresponding extremal co-dimension two surface in the 
bulk. In the following, without loss of generality, we set Newton’s gravita¬ 
tional constant Gn = 1 . 

These expressions are applicable to all holographic models. For AdS 4 , 
which is of interest here, dA is a closed curve and H extr is two-dimensional. 
Then, the area of the extremal surface, which itself will be denoted by H extr 
in the following, is expressed in terms of the induced metric as 



A extr = [ d 2 ay/ 7 , 

( 2 . 8 ) 

where 

dX» (a) dX v (a) 

lab ~ 9 au d(ja d(jh 

(2.9) 

and 


7 = det (lab) ■ 

( 2 . 10 ) 


Here, g is the bulk metric, a a are coordinates parametrizing the extremal 
surface, X (a) are the parametric equations of the extremal surface in the 
bulk and 7 is the induced metric on the extremal surface. 
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2.3. Entanglement Entropy in Global AdS . 4 for a Polar Cap Region 

Extremal surfaces in AdS space-times have been mostly studied in Poincare 
patch coordinates. I11 those coordinates, the space-time line element of AdS 4 
with unit scale is 

ds 2 = (— dr 2 + dz 2 + dx 2 + dy 2 ) . (2.11) 

z 

Then, for a disc region A in the boundary plane described by (x — Xq )" + 

(y — y 0 ) 2 < R 2 , the corresponding extremal surface A extr in the AdS bulk is 
given by 


(x - x 0 f + (y- y 0 ) 2 + £ 2 = R 2 , 

r = r 0 . 


( 2 . 12 ) 


Without loss of generality we may take the disc centered at (0, 0). 

Passing from Poincare coordinates (r, z, x, y) to global coordinates (t, r, 9, f) 
with the aid of the coordinate transformation 


r = 


\Jr 2 + 1 sin t 


(2.13) 


\/r 2 + 1 cos t + r cos 9 ’ 

1 

\Jr 2 + 1 cos t + r cos 6 
r sin 9 cos ip 

Vr 2 + 1 cos t + r cos 9 ’ 
r sin 9 sin ip 

\Jr 2 + 1 cos t + r cos 9 
the space-time metric takes the form 

ds 2 = - (r 2 + l) dt 2 + j- + r 2 (d^ 2 + sin 2 6>d^ 2 ) , (2.14) 


z = 


x = 


V = 


whereas the extremal surface (2.12) corresponding to the choice (a; 0 , Vo) = 
(0, 0) is given in global coordinates by 


t = t 0 , r ( 9 ) = 


cos 0\/tan 2 0o — tan ^9 


9 e [0, 6»o], <pe [0, 2tt) , (2.15) 
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where 9 0 and t 0 are specific functions of R and r 0 
details). Equivalently, we have parametrization 


see 


Appendix A for more 


t — t 0 , 


9 (r) — arccos 



r G [cot 9 0 , oo), (p G [0, 2tt) . 


( 2 . 1 %) 

The region A becomes a polar cap in global coordinates and the complemen¬ 
tary cap is the region A c . 

Figure 1 depicts the regions A and A c on the spherical boundary of space- 
time. 



region A c 


Figure 1: The region A has opening angle 9q. The entangling curve dA is the thick line. 


It is easy to confirm that the surface (2.16) is extremal, obeying the 


particular conditions derived by minimizing the area functional (2.8), under 


the assumption that the surface is rotationally symmetric, i.e., t — t (r) and 
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9 = 9 (r), since 


(r 2 + 1) t'r sin# 


(r 2 + 1 ) f ' 2 + + r 2 #' 2 

r 3 0' sin 9 

(r 2 + 1) f' 2 + + r 2 9' 2 


= 0, 


= rcos9\ — (r 2 + l)f /2 + 


(2.17) 


r 2 + 1 


+ r 2 9 ' 2 , 


(2.18) 


where prime denotes differentiation with respect to r. 

The extremal surface emanates from a polar cap boundary region de¬ 
scribed by 9 < 9q. An introduction of non-vanishing parameters Xo or yo 
would rotate the entangling curve so that its symmetry axis would not any¬ 
more coincide with the axis corresponding to the azimuthal angle q h Having 
said that, we stick to the choice xq = yo = 0 from now on. 

Figure 2 depicts the extremal surface A extr that emanates from a polar 
cap region in the boundary of space-time in global coordinates and extends 
in the interior of space-time. The radial coordinate in the plot is proportional 
to the so called tortoise coordinate, arctanr. 



Figure 2: The extremal surfaces in AdS^ space-time for various choices of 9q. 






















Parametrizing the extremal surface with the coordinates 0 and r, the 
induced metric turns out to be diagonal with elements 


7rr 


sin 1 2 * 0 n 


r 2 (l + X) [l - (l + cos 2 0 o ] ’ 
1 


= r 


1 + 0T ) COS 2 6»o 


(2.19) 

( 2 . 20 ) 


The determinant of the induced metric is 7 = sin 2 0 o /(l + l/'r 2 ) and the area 


follows from equation ( 2 . 8 ), 


y4 extr = 2 tt lim (r sin 0 O - 1 ) 


( 2 . 21 ) 


The hrst term is the divergent “area law” term, while the second one is 
universal independent of the UV cutoff [30, f3Tj . 


2-4- Perturbations of Entanglement Entropy 

The area functional (2.8) depends on the background metric g as well as 
on the embedding variables A", and, of course, it also depends implicitly on 
the entangling curve dA curving the region A in the boundary held theory. 
Perturbations of the bulk metric induce changes of the minimal surface. The 
variation of its area is given in general by 


M extr 


(L4 extr (g, X) 
Tg 


Sg + 

9=90, X=Xq 


<h4 extr (g, X ) 

Jx 


5X 


9=90, X=Xq 


(2.23) 


as both the metric and the embedding equations of the extremal surface vary 
around their unperturbed values go and A 0 . 


1 Recall that the entanglement entropy for a disk region of radius R in the boundary 
of AdSd+i is given by 


ai{R/L) -\-aAR/L ) + ■ • • + 2 {P/E) + a d-i, d odd, 

ai(R/L) d_2 + a 3 (R/L ) d ~ 4 + ... + a d - 3 {R/L f + a d -dog (R/L ), d even, 


( 2 . 22 ) 

where L is the UV cutoff. The first term is the “area law” term. For d even, the logarithmic 
term is universal and connected to the conformal anomaly 0 HDl M IM32]. For d odd, 
which is the case of interest here (d = 3), the constant term is universal and it obeys a 
holographic “c-theorem” mm- 
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When the perturbation obeys Dirichlet boundary conditions, the region 
A and its boundary dA remain fixed. Then, the second term vanishes, as the 
original surface described by equations Xq extre mi ties the area functional with 
the given loop dA. When the metric perturbations satisfy general boundary 
conditions, computing the variation of the area of the minimal surface is 
more subtle. For this, we suppose that the boundary loop dA specifying the 
extremal surface is determined by a set of parameters b, which may also vary 
as metric perturbations are turned on. Then, the area functional has the 
particular form 

^extr = ^extr ^ x ^ ^ (2.24) 

and its variation breaks down as follows, 


5A extr 


S yfextr X ) 


Sg 


Sg 


+ 


9=90, X=Xo 

5A extl (g, X ) 


SX 


SX(g,b) 


9=90, X=X 0 

5A extI (, g , X) 


^g 


Sg 


sx 


9=90, b=bo 

dX (g, b ) 


9=90, X=X 0 


db 


5b. 


9=90 ,b=b 0 


(2.25) 


The second term vanishes, as it corresponds to the variation of the area for 
variations of the extremal surface with a fixed boundary that is provided by 
the entangling loop dA. The last term can be simplified if one considers the 
area as functional of the metric as well as the parameters b specifying the 
entangling loop. As a result, the total variation of the area takes the form 


5 A = 


(L4 extr (g, X) 


Sg 


Sg + 


dA extr (g, b ) 


:= 6A extr + 


9=90, X=Xo 
extr 


db 


5b 


9=90, b=bo 


(2.26) 


hereby defining the individual contributions 5A extr and 5Al xtr . 

The first term follows by varying the relation (2.8) and the result is written 
in terms of the unperturbed induced metric y 0 and 5~f a (, as 


5A 


extr 

9 


1 

2 


/ 


d 2 cr \/7o(7o) aft ^7a6, 


(2.27) 
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where 


(7o ) ab ~ (do) 
hob = Sg lw 


dXf {a) dX Q v (a) 

da a da b ' 
dX 0 » {a) dX 0 " (a) 


da° 


da b 


(2.28) 

(2.29) 


The second term can be easily calculated for polar caps, which only depends 
on a single parameter b that is taken to be 9 0 . Then, SAl xtT follows by varying 


the relation (2.21) with respect to the cap parameter 9 0 . 


The individual terms 5A extr and 5Al xtI will be explicitly computed later 
for small perturbations of the AdS^ metric satisfying the linearized Einstein 
equations with general boundary conditions. 

2.5. Modular Hamiltonian for a Polar Cap Region 

Unlike the Ryu-Takayanagi formula for expressing the entanglement en¬ 
tropy in terms of the bulk gravitational theory, there is no similar expression 
for the modular Hamiltonian. In general, the modular Hamiltonian is a non¬ 
local operator and there is no way to this day to find an expression for the 
modular Hamiltonian for a general boundary state and region A. In some 
cases, however, the modular Hamiltonian generates a geometric flow leading 
to a local expression. 

Such an example is provided by taking a disk region A of radius R in 
a Minkowski space boundary [IS] , in which case the modular flow in the 
Cauchy development of the disk can be connected with the modular flow in 
Rindler space through a conformal transformation; the details are provided 
Appendix B| Then, the modular Hamiltonian is expressed in terms of the 


m 


holographic energy-momentum tensor and a conformal Killing vector that 
leaves invariant the entangling curve dA and its causal development, as 


SE a = [ dWT^C, (2-30) 

Jc 

where C is a space-like surface with boundary dA and dZU is the differential 
volume form on C HZ]. If C is selected to be a constant time slice, it will 
coincide with the region A. For a disk in Minkowski space, this conformal 
Killing vector ( is a linear combination of the Killing vectors corresponding 
to the time translation and a special conformal transformation; it is selected 
so that ( vanishes on dA and its causal development. 
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Using the coordinate transformation (2.13), the conformal Killing vector 
in question takes the following form in global coordinates (for more details 
see 


Appendix A), 


c = 


2?r 
sin 6 r 


[(cos ( t — t 0 ) cos 9 — cos 9 0 ) d t — sin (t — t 0 ) sin 9d e ] , (2.31) 


which indeed satisfies the conformal Killing vector equation 


+ V,C M - ~ d 9^xC, X = 0 (2.32) 

and at the same time vanishes at the entangling surface 6 = 9 0 and t = t 0 - 
The conformal Killing vector ( is the boundary limit of a conformal Killing 
vector held £ in the bulk, 


£ = 


2tt 

sin 9n 


r cos 6 cos (t — t 0 ) 
a Jr 2 + 1 


- cos 6* 0 d t 


/ 2 , Q ■ U + \ n Vr 2 + lsm9sm(t -t 0 ) 

+V r 2 + 1 cos 9 sm ( t — t 0 ) o r - Og 


, (2.33) 


which also vanishes at the entangling surface. 

Note in passing that if one considers generalized theories of gravity (i.e., 
higher derivative corrections) an appropriate generalization of Ryu-Takayanagi 
conjecture will be required. Such generalizations can be obtained through 
the Iyer-Wald theorem that involves bifurcate Killing horizons generated by 
the conformal Killing vector £ given by (2.33) [20]. 

Formula (2.30) can be evaluated for any space-like surface with boundary 
identical to the entangling surface. Selecting the t = t 0 surface and using the 
specific form of the Killing vector (2.31), the modular Hamiltonian for the 
polar cap region 9 < 9q is written as 


4tt 2 f 9 ° 

8Ea —— —t p / d6*sin6 l (cos6* — cos0 o ) 7(0 (2.34) 

sm 6*o Jo 


This completes the presentation of the general formulae for the quantities 
of interest. In section [4] 5E will be calculated and compared to 8 S for all 
different kinds of AdS^ perturbations satisfying general boundary conditions. 
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3. Linearized Gravity in AdS 4 

Now we come to the theory of gravitational perturbations of AdS^ space- 
time, allowing for general boundary conditions. It is convenient to split the 
perturbations into two complementary sets with opposite parity and reduce 
the linearized Einstein equations to an affective Schrodinger problem, which 
is exactly solvable. Then, we write down the holographic energy-momentum 
tensor for all such perturbations, which will be used later. The material we 
present here is based on the discussion found in [26] (but see also [27] for an 
overview of the subject) and references therein. 


3.1. Linear AdS 4 Perturbations 

We consider Einstein gravity in four space-time dimensions with negative 
cosmological constant A, 

R^ = Ag^. (3.1) 

AdSi is the maximally symmetric solution whose metric takes the following 
form in spherical coordinates, 

n r 2 

ds 2 = —f (r) d.t 2 + + r 2 idO 2 + sin 2 #^ 2 ) , (3.2) 

/ ( r ) v 


where 


/ M = i - If 2 . 


(3.3) 


We set the AdS scale a/— 3/A equal to one to simplify the presentation. 
For later use, we define the tortoise coordinate x as 


dx = 


dr 

fir) 1 


which is given explicitly as 


r = tan x . 


(3.4) 


(3.5) 


The tortoise coordinate is an angular variable ranging from 0 to 7r/2 as r 
varies from 0 to oo. 

Next, we consider linear perturbations around the AdS 4 metric (3.2) satis¬ 
fying Einstein’s equations. The gravitational perturbations fall in two com¬ 
plementary classes with opposite parity, namely axial and polar perturba¬ 
tions. Without loss of generality we restrict attention to axially symmetric 
deformations, which are described in terms of Legendre polynomials Pi (cos 6 ) 
instead of more general spherical harmonics Yj m ( 6 ,<p). 
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Axial Perturbations 

The metric perturbations of this class are parametrized by two functions 
ho (r) and hi (r) as follows, 


ds 2 = —/ (r) dt 2 + 


dr 2 

f(r) 


+ r 2 (< id 2 + sin 2 dd</? 2 ) 


2e lut sin 6 —' " —- (d 0 (r) dt + /q (r) dr) d<^, (3.6) 

dO 


up to reparametrizations. It turns out that Einstein’s equations are equiva¬ 
lent to the effective Schrodinger problem with respect to the tortoise coordi¬ 
nate, 


d 2 T a 


(x, 


dx 2 




sin 2 a: 


(x) = orT* 


(x 


(3.7) 


where the functions ho (r) and hi (r) are expressed in terms of the solutions 
as 


ho (x) 
hi (x) 


id . .. 

(tanx W ax (x)), 

co dx 

sin x cos x T ax (x). 


(3.8) 

(3.9) 


Polar Perturbations 

The metric perturbations of this class are parametrized by three functions 
H 0 (r), Hi (r) and K (r) as follows, 


ds 2 = —f (r) dt 2 + f r + r 2 (d,9 2 + sin 2 9dtp 2 ) 

J (r) 


+ e~ lwt Pi (cos 9) 


Ho (r) ( / (r) dt 2 + 


dr 2 

f(r) 


2 Hi (r) dtdr + K (r) r 2 (dd 2 + sin 2 dd99 2 ) 


, (3.10) 


up to reparametrizations. Similarly to the previous case, Einstein’s equations 
for polar perturbations turn out to be equivalent to the effective Schrodinger 
problem with respect to the tortoise coordinate, 


d 2 T 


pol 


[X, 


dx 2 


+ d!±iT po i 


sin 2 a; 


x 


= orT 


pol 


[X 


(3.11) 
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which is identical to the one describing the axial perturbations. The functions 
Hq (r), H i (r) and K (r) are expressed in terms of the solutions of the effective 
Schrodinger problem as 


H 0 (x) 
H^x) 
K (x) 


l (l + 1 ) 2 ■ 2 d \ t / \ 

-cot X — UJ Sin X COS X + cos x — Wpol ( X) , 

2 ax. 

—iu cos x— (sinx^/p 0 i (x)), 
ax 


1(1 + 1 ) 


cot x + -^j Tpol (x). 


(3.12) 

(3.13) 

(3.14) 


Thus, the spectrum and the eigen-functions of the operator —d 2 /dx 2 + 
1(1 + l)/sin 2 x in the closed interval x G [0, n/2] completely determine the 
gravitational perturbations of AdS 4 space-time. The general solution of the 
linearized Einstein equations is written as linear combination of the axial and 
polar solutions. 


3.2. Boundary Conditions 

The solution of effective Schrodinger problem can be expressed in terms 
of hypergeometric functions. The normalizable solution, which vanishes at 
r = 0, is 


T (x) = cosxsin 1+1 x F ( - (l + 2 + u >), - (l -f 2 — w); l + sin 2 x ) . 


(3.15) 


The other independent solution of the effective Schrodinger problem is 


cosx 


sm ; x 


= — r F ni-t + t+u), -(~l + l-u);^ -Z; sin x , (3.16) 


but it diverges as r —> 0 and is not normalizable. 

The boundary conditions as r —> oo can be systematically described by 
expanding the normalizable solutions in powers of 1/r, as 




/o + 



r 




(3.17) 


Using properties of the hypergeometric functions, it turns out that I 0 and I\ 
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are expressed in terms of Gamma-functions as 


r- 1 (i(/ + 2-^)) , (3.18) 

h — -2T 1 f-(/ + 1 + w)i T 1 f-(/ + 1 - ~\i) . (3.19) 

The coefficients of all other terms are expressed in terms of Jo and /1, but 
the details do not really matter. We only note here, for later use, that 

J 2 = |(/(/ + l)-u; 2 ). (3.20) 

The boundary conditions imposed as r —> oo are solely described in terms 
Jo and ii, since 


/n = 


r 1 ( — (/ + 2 + lS) 


J 0 = T(r = oo).* h = - — (r = oo). (3.21) 

Thus, solutions, characterized by Jo = 0 obey Dirichlet boundary condi¬ 
tions for the effective Schrodinger problem, while solutions characterized by 
Ii = 0 obey Neumann boundary conditions. Solutions with Jo and I\ taking 
more general values correspond to more general (mixed) boundary conditions 
determined by the ratio Iq/I\. In the following, we will denote by 1^ the co¬ 
efficients of the large r expansion of T(r) associated to axial perturbations 
and Jk for the polar perturbations, for distinction. 

For given boundary conditions, the spectrum is discrete; for example, for 
Dirichlet boundary conditions we find 

u) D = ± (2 n + 1 + 2), (3.22) 

whereas for Neumann boundary conditions the spectrum is 

co N = ± (2 n + 1 + 1). (3.23) 

More general boundary conditions give rise to intermediate (generally not 
equidistant) frequencies. In all cases, the spectrum is implicitly determined 
by ratios of Gamma-functions fixed by the ratio Iq/Ii or Jq/J\. 
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3.3. The Holographic Energy-Momentum Tensor 

The energy-momentum tensor is divergent in asymptotic AdS spaces and 
appropriate renormalization is required to make sense of it, using counter¬ 
terms. Here, we summarize the results of holographic renormalization [33ff36| 
for the boundary space-time metric and the corresponding energy-momentum 
tensor for gravitational perturbations of AdS^ space-time satisfying general 
boundary conditions, following |26j . Thus, we have per sector the following 
results, setting Newton’s constant Gn = 1: 

Axial Perturbations 

The boundary metric after conformal rescaling takes the form 


ds 2 = — dt 2 + (dO 2 + sin 2 0dip 2 ') + luJt sin 6 ^ ^ dtdip , (3.24) 

v ’ co do 


whereas the non-vanishing components of the energy momentum tensor for 
small perturbations are 



(3.25) 


(3.26) 


Polar Perturbations 

In this sector, the boundary metric after conformal rescaling takes the 
form 


ds 2 = — d,t 2 + (l — Jie luJt P l (cos 9)) (d6 2 + sin 2 6 d(p 2 ) , (3.27) 
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whereas the non-vanishing components of the energy-momentum tensor for 
the corresponding small perturbations are 


8 v iTu = J (Z - 1) l (l + 1) (l + 2 ) Joe-^Pi (cos 9) , (3.28) 

8 t rTee =-h {l + 1) (l - a; 2 ) J 0 e~ iujt Pi (cos 9) 

- j (l (l + 1 ) - 2 w 2 ) J 0 e-^ cot o dPl (™ s9 ) ? ( 3 . 29 ) 

4 at) 

8nT w = (Z + 1 ) (Z (Z + 1 ) — 1 — cu 2 ) Joe _ ^sin 2 0 Pi (cos 6 *) 

+ ^ (Z (Z + 1) — 2cu 2 ) J 0 e~ lut sin 9 cos 9 ^ 1 , (3.30) 

8tt7^ = (Z - 1) (Z + 2 ) J &~™ — . (3.31) 


In all cases, the energy-momentum tensor is traceless and covariantly 
conserved, 

T a a = 0 , V a T ab = 0 (3.32) 

with respect to the boundary metric, as required on general grounds. 


Equations (3.24) and (3.27) imply that for axial perturbations Dirich- 


let boundary conditions for the effective Schrodinger problem correspond 
to Dirichlet conditions for the metric, while for polar perturbations Neu¬ 
mann boundary conditions for the effective Schrodinger problem correspond 
to Dirichlet conditions for the metric. These results will be particularly useful 
while computing the modular Hamiltonian as an appropriately weighted inte¬ 
gral of the energy density for all kind of perturbations. It can already be seen 
that SE always vanishes for axial perturbations, since Ttt = 0 irrespective of 
boundary conditions, whereas SE does not vanish for polar perturbations 
satisfying boundary conditions other than Jq = 0 . 


4. Entanglement Entropy for Gravitational Perturbations 

We are now in position to compute holographically the variation of the 
entanglement entropy and the modular Hamiltonian for gravitational per¬ 
turbations of AdS<x space-time satisfying general boundary conditions. First, 
we compute the variation of the area of the extremal surface emanating from 
the loop dA , as 5H extr = 8A e * tT + 5H^ xtr , then we compute SE, and, finally, 
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we compare the results with the first law of thermodynamics for the entan¬ 
glement entropy, 5S = 5E. We find that the deformations of the entangling 
curve dA should be isoperimetric, when the boundary conditions of the met¬ 
ric allow for fluctuating geometries on the boundary of AdS^ for, otherwise, 
the first law of thermodynamics would not be realized holographically. Our 
results generalize earlier work on the subject, going from Dirichlet to more 
general boundary conditions for the metric. 

4-1. Calculation of 5A e * tr 

First, we work out the contribution 5A e * tr to the variation of the area 
of the minimal surface. Naturally, we consider the effect of axial and polar 
perturbations separately. 


Axial Perturbations 

For axial perturbations, the variation of the induced metric on the mini¬ 
mal surface, which is generally given by equation (2.29), turns out to be 

(4.1) 


S^y rr 0 5 


d 


hrp = $g r <p = h x (r) e luJt ° sin#— P L (cos 9) 


6 =arccos (cos 


5 0O\/1+4t) 


(4.2) 


and so the corresponding variation of the area, given by equation (2.27), is 

SA e * tT = 0 . (4.3) 

Polar Perturbations 


For polar perturbations, the variation of the induced metric (2.29) turns 
out to be 

c . fdd\\ 

O^frr — og rr -f- I J dggg 


1 + r 2 


Ho (r) + 


cos 2 0 o A' (r) 


r 2 (1 - (1 + y cos 2 #o) 


e-<*iqcos<wi + y ], 

(4.4) 


y, = Sg vv = r‘ 1 - 1 + -r cos 2 0 o I< (r) e cos0 o \ 1 + -5 , 




(4.5) 

(4.6) 
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Then, using equations (2.19), (2.20) and (2.27), we find after some algebra 


that 


5A e ™ = 7T 


r ^ 

cote 0 ^/(l + T) 


(Hq (r) + K (r )) 


r' 2 tan 2 0 o 


(K (r) — H 0 (r)) 


e-^P, cos0o\/l + ^ ■ (4-7) 


Trading r with the tortoise coordinate x and expressing Hq and K in terms 


of the Tp 0 i, as given by equations (3.12) and (3.14), we obtain 


SAf* = 7T 


dx 


sin d 0 


1 - 0o 


cosx 


+ tan x 1 + 


l (l + 1) d'poi — y — 
cos 2 a; \ dTpoi 


sin 2 dr 


cos 2 x 

sin 2 d 0 

—iujto 


co 2 ^ 


dx 




pol 

COS do 
sin x 


. (4.8) 


The effective Schrodinger problem (3.11) can be used further to eliminate u ;, 
leading to the following expression, 


8Af* = 7T 


/* 

" d 

f sin d 0 cos x \ 

dty pol 

dx 

\ cos x sin d 0 ) 

dx 



+1 (Z + l) 


COS do cot X 




pol 


-ILjto p 


COS dr 


Sill X 


(4.9) 


tan d 0 sin x 

ft is convenient at this point to employ the following identities among Leg¬ 
endre polynomials, 


4 (P m (f) - fS_, «)) = (2 1 + 1) P, (£), 


(4.10) 


4p, K) = (p, + , K) - p,_, «)). (4.11) 


to perform integration by parts and get 
"sin d 0 


8A e p = ire~ iut0 


cos x sin dn 


cos x \ rfT po i f cos do 

Pi 


dx 


sin x 


J f -00 


— 7re 




cot 9 o hi±lifp 1+ / cos,, » 


2Z + 1 


sm a; 


-Pi-1 


COS dr 


sm a; 




pol 


f-0o 

(4.12) 
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The contribution from the end-point x = | — 9 0 vanishes for both bracketed 
terms, and, thus, the result is 


6A extl = Tre~ iujt ° lim 
y ,7r 


cot 9 1 


sin^o dTpoi 
cos x dx 
1(1 + 1 ) 


Pi (cos 0 O ) 


(Pi+ 1 (cos 6*o ) - P /-1 (cosdo)) V I ; pol 


(4.13) 


21 + 1 

Going back to the coordinate r and using the asymptotic expansion of the 
Schrodinger wave-functions (3.17) together with equation (3.20), we arrive 
at the final expression 


§A extr = _ ne -iut 0 


lim (rJi + J 0 (l (Z —(— 1) — to 2 )) sin d 0 Pi (cos d 0 ) 

—>oo v v '' 

1(1 + 1 ) 


+ cot 9 0 


21 + 1 


-J 0 (P t+ 1 (cos 9q) - Pi- 1 (cos^q)) 


Note that 5A e * tr is divergent when the perturbations of the metric 
satisfy Dirichlet boundary conditions, J\ ^ 0. 


(4.14) 
do not 


4-2. Calculation of 5Al xtT 

Since we are considering axially symmetric perturbations of the bulk met¬ 
ric, the deformations of the entangling curve are taken to be independent of 
the angular coordinate <f> to conform with the symmetry. Thus, the variations 
of the parameter 9 0 in the polar cap are taken to be independent of 0, i.e., 


(t) — 9 q + 69 ( t ). 


(4.15) 


Consequently, 6A^ tr can be derived directly by varying (2.21) with respect 
to do, leading to the following result 


6A c b x r = 27 t lim r cos 9 0 69 (t 0 ) 


(4.16) 


for all kind of metric perturbations. 

Note that 6A^ tr also appears to be divergent on general grounds. It 
is important to see how these divergences combine with those appearing in 


6A e * tr to form 5A extr . 
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4-3. Net Variation of the Entanglement Entropy 

Summing up the two contributions, we obtain the net variation of the 
area of the minimal surface, <5A extr , and, hence, the variation of entanglement 
entropy, 5Sa • We present the result for axial and polar perturbations, 

setting Newton’s constant Gn = 1. Thus, for axial perturbations, we have 

^axiai = lim r cos do 50 axial ^ , ( 4 . 17 ) 

2 r—»oo 

whereas for polar perturbations the result reads 


5^ olar 


- cot 0 O + r ,7 ° ^ Pl+l ( cos6, °) 


to) - n- 1 fCOS C'oJ J 
J lim rcos0 o 50 polar (to). (4.18) 

2 r —>oo 


The variance of the entropy appears to be divergent, which is unphysical. 
Specific conditions should be imposed on <50 polar (to) to make them cancel, but 
these will be discussed later, while comparing SSa with SEa, together with 
their geometric meaning. 


4-4- Modular Hamiltonian for Gravitational Perturbations 

Next, we compute the variation of the modular Hamiltonian by specializ¬ 
ing formula (2.34) to axial and polar perturbations satisfying general bound¬ 
ary conditions. As before, Newton’s constant is set equal to 1. 


Axial Perturbations 

Since E t t — 0 for the axial perturbations, we immediately obtain 

8 Ef VAl = 0 (4.19) 

in all cases, irrespective of boundary conditions. 
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Polar Perturbations 

For polar perturbations with Jo 7 ^ 0, the component 7tt of the holographic 
energy-momentum tensor does not vanish, and, hence, the corresponding 
variance of the modular Hamiltonian differs from zero. 

First, substituting the expression (3.28) for % t into equation (2.34), yields 
the following expression for the variance of the modular Hamiltonian, 

r -00 


= ir(J-l)Ui+l)(i + 2) 

8 sin 6 *o 


dd sin 9 (cos 9 — cos 6 * 0 ) Pi (cos 9). 

(4.20) 


Then, we perform integration by parts and make use of formula (4.10) twice 
to find 

££ P olar = 7T (/-!)/(/ + !)(/ + 2) 


8 sin 6* 0 

x 


21 + 1 


Pi+ 2 (cos 6 *q) Pi —2 (cos 6*o) 2 (21 + 1) P t (cos 6* 0 ) 


21 + 3 ' 21-1 (21 - 1) (21 + 3) 

Next, using the following identity of Legendre polynomials, 

(l + 1) Pi +1 (0 + l Pi— i (0 = (2/ + 1) £ P, (0 

to eliminate P ;_2 and -Pz+ 2 , we arrive at the expression 


(4.21) 


(4.22) 


rppolar _ 7T (/ — 1) l (l + 1) (l + 2) 
" “ ~ ' " 21 + 1 


X 


8 sin 6* 0 
fP 

cos 6*0 


Joe 


— 


P/+i (cos 6*o) P /-1 (cos 6*o) 


1+2 l-l 

which can be conveniently rewritten as 


(21 + 1 ) Pi (cos 6 *q) 

(l - 1) (l + 2 ) 


(4.23) 


r p polar _ 7T (/ - 1) l (l + 1) (l + 2) 

- 8ita^ 2lTl J ° e 

((i + 1) P;+l (cos 6 *q) + ZP|_i (cos^n)) 


X 


L(/ - 1 ) (/ + 2 ) 

2 cos 6*0 


(Z -!)(/ + 2 ) 


(-Pz+i (cos 6 * 0 ) - Pi—i (cos 6 * 0 )) 


( 2 / + 1 ) Pf (cos 6 * 0 ) 

(l - 1) (l + 2 ) 


(4.24) 
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Using once more the identity (4.22) yields the final result, 


5E P ° lar = —— e ~ iu}t ° 


-Jo 1(1 + 1) sin 6 0 Pi (cos 0 O ) 


+ cot 0 O + Jq (Pl +1 (cos 00 ) - Pi -1 (cos 0 O )) 


. (4.25) 


2Z + 1 

We see that in all cases the variance of the modular Hamiltonian is finite. 


4-5. Comparison with First Law of Thermodynamics 

The first law of thermodynamics SS = 5E for the entanglement entropy 
and the modular Hamiltonian is a tautology from the point of view of the 
boundary theory. Its realization through the Ryu-Takayanagi formula (2.7) 
and the holographic energy-momentum tensor is not automatic, but it is a 
consistency check of the whole scheme. Comparing the results derived above, 
we find that 5Sa = 5Ea for any polar cap region A, provided that the en¬ 
tangling curve at the boundary also deforms, where appropriate, in a certain 
way. The boundary is not held fixed when general boundary conditions are 
imposed on the metric, and, thus, the polar cap region A and the entangling 
curve dA also vary. Thus, 59(to) should be adjusted accordingly to achieve 
agreement with the first law of thermodynamics. As will be seen shortly 
this requirement has natural geometric interpretation. It also gets rid of the 
divergences appearing in the computation of SSa- 


Axial Perturbations 

First, we consider the case of axial perturbations and compare the ex¬ 
pressions (4.17) for SSa and (4.19) for 8 Ea- The first law is realized holo- 

<50 axial (t 0 ) = 0 (4.26) 


graphically provided that 


irrespective of boundary conditions. Thus, axial perturbations do not require 
adjustment of the polar cap region. 

A way to understand this is to consider the boundary metric of the per¬ 
turbed AdS/i space-time, prior to conformal scaling, and work out the two- 
dimensional induced metric on the corresponding polar cap region. For axial 
perturbations it reads 


<^s 2 axia i = r 2 (d9 2 + sin 2 9d(p 2 ') . 


(4.27) 


The polar cap is placed on a round sphere that remains inert to perturbations 
irrespective of boundary conditions. 
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Polar Perturbations 

Next, we consider the case of polar perturbations and compare the cor¬ 
responding expressions (4.18) for 5 Sa and (4.25) for SEa ■ The first law of 


thermodynamics for entanglement entropy is realized holographically pro¬ 
vided that the entangling curve in the polar cap region varies according to 
the equation 


<56 |polar (t 0 ) = -e~ iut ° 


r \ 2 


tan 9 0 Pi (cos 9 0 ). (4.28) 


This also takes care of the unphysical divergences appearing in 5Sa■ Note at 
this end that a subleading term proportional to J 0 is included in the variation 
59 polaT (t 0 ); it contributes a finite term to SSa which should also be removed 
to match it with SEa- 

In this case, perturbations of the bulk metric require adjusting the en¬ 
tangling curve in the boundary. To understand what is happening, let us 
consider the boundary metric of the perturbed AdS± space-time, prior to 
conformal scaling, and work out the two-dimensional induced metric on the 
corresponding polar cap region A. For polar perturbations, the induced met¬ 
ric on the spherical slices of the boundary is 


ds 2 poiar = r 2 (l + K (r) e lut Pi (cos 9)) ( d9 2 + sin 2 9d(f 2 ) , 


(4.29) 


where Ii(r ) is the radial function appearing in the general parametrization 


(3.10) of the bulk space metric. The asymptotic form for K (r) follows from 
the large r expansion of the Schrodinger wave function T po i, via equation 


(3.14), and it reads 


Jo 


K (r) = - Jr - — -Z (l + 1) - u; 2 + O ^ 


(4.30) 


It implies, in particular, that the condition (4.28) on 59 polar (t Q ) takes the 
following equivalent form 


39 polar (t 0 ) = — -e~ iujt °K(r) tan 9 0 Pi (cos 9 0 ) 


(4.31) 


up to irrelevant subleading terms. 

The deformation of the entangling curve is such that its line element is 


inert to the perturbation. Using the spherical metric (4.29), we find that 
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the induced line element of the boundary curve of a polar cap region with 
parameter 9 is 


dl 


lim J~g^dip 

r—>■ oo 


lim 

r—> oo 


! + - Jt (r) e 


—iujt ; 


.P; (cos #) I r sin 9 dip . 


(4.32) 


Here, we also take the large r limit as the computations take place on the 
boundary of space-time, prior to rescaling. Setting 6 = 9 0 + 59(t 0 ), where 
69(to) is provided by (4.31), we have to linear order that 

(l + \K (r) e~ lut °Pi (cos (9 0 + 59 (f 0 )))) r sin (9 0 + 59 (t 0 )) = 

(l + \K (r) e~ zuto Pi (cos 6 * 0 )) r sin (9 0 + 59 (t 0 )) = rsm9 0 + O (I) (4.33) 


Thus, although the polar cap region and its boundary undergo perturbations, 
the line element of the entangling curve remains constant, 


dl = lim r sin 9q dip . 

r—>• 00 


(4.34) 


Summarizing, we see that for all kind of gravitational perturbations of 
AdS 4 space-time (axial or polar) and for all kind of boundary conditions, 
the holographic realization of the first law of thermodynamics requires that 
the line element of the entangling curve remains invariant. Region A, whose 


entanglement entropy is provided by formula (2.7), deforms with time in 


such a way that the line element of the entangling curve remains constant. 
This requirement also takes care of the divergencies that otherwise would 
have inflicted the variance of the entanglement entropy. Since the first law 
of thermodynamics is trivially valid from the CFT point of view, the above 
statement should be viewed as addendum to Ryu-Takayanagi prescription. 

Although the laws of entanglement thermodynamics are not yet com¬ 
pletely understood, it has been shown that for spherical regions the reduced 
density matrix pa is identical to a thermal density matrix with effective in¬ 
verse temperature equal to the perimeter of the entangling circular curve 
[ IS] . Through this connection, our results can be thought to provide the 
holographic interpretation of an isothermal process of entanglement thermo¬ 
dynamics. 

An interesting question arises by comparing our calculation to previous 
ones [19) 20] - In those papers, the first law of thermodynamics for entangle¬ 
ment was realized holographically using gravitational perturbations satisfy¬ 
ing Dirichlet boundary conditions, without making reference to the contribu¬ 
tion 5Al xtl . In our case, even for Dirichlet boundary conditions for the polar 
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perturbations of the metric, J\ = 0, 8Al xtl is non-vanishing. Its contribution 
is absolutely necessary for validating the first law of thermodynamics for en¬ 
tanglement. The resolution to this ostensible contradiction is that each one 
of the terms arising in 5A extT and 8A e ^ tr depend on the choice of coordinates, 
but their net sum is, in fact, coordinate independent when the line element 
of the entangling curve remains constant. The variation <5A^ xtr vanishes for 
Dirichlet boundary conditions only in some coordinate systems, which in¬ 
clude the Fefferman-Graham coordinates used in references mm- More 
technical details about this issue are provided in |Appendix C 

The demand for isometric deformations of the entangling curve is also 
strongly connected with renormalization issues of the entanglement entropy. 
In odd dimensional conformal field theories, as in our case, the finite contri¬ 
bution to the entanglement entropy can be contaminated by the divergent 
terms, using alternative definitions of the UV cutoff scale L. This issue has 
already been noticed in the literature [T8J. The universal constant term can 
be isolated using an appropriate renormalized entanglement entropy, which 
in 2 + 1 conformal field theories (and at least for special choices of the en¬ 
tangling curve) takes the form 133 

See :=S EE - J^- L{dA), (4.35) 


where L ( dA ) is the length of the entangling curve, resolving the above is¬ 
sue. This is also in agreement, via the Ryu-Takayanagi conjecture, with the 
definition of renormalized area in the bulk theory, provided by Graham and 
Witten [38] : in that work the same contamination arises from the divergent 
terms, resulting in a specific choice of radial coordinate to cure the prob¬ 
lem. When the metric perturbations of AdS^ are taken in the form (3.6) or 
(3.10), the corresponding radial coordinate r does not belong to that special 
class, calling for a renormalized definition of holographic entanglement en¬ 
tropy, as in equation (4.35). In this language, the term 8A extr corresponds to 
the variation of See , whereas the term 8Al xtI corresponds to the variation 
of the counter-term. Then, as it is implied by our results, the inclusion of 
the counter-term is equivalent to the statement that the unrenormalized en¬ 
tanglement entropy, which is provided by the Ryu-Takayanagi formula (2.7), 
is computed for a region A bounded by an entangling curve that deforms 
isometrically. 
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5. Electric-Magnetic Duality for Entanglement Entropy 


In this section, we review briefly the electric-magnetic duality exhibited 
by linearized gravity, together with its holographic manifestation, following 
[26], and examine the implications to entanglement entropy. 


5.1. Energy-Momentum/Cotton Tensor Duality 

Linearized gravity around maximally symmetric backgrounds exhibits a 
rank-2 generalization of electric-magnetic duality in four space-time dimen¬ 
sions [221 [23]. It is best described in terms of the Weyl tensor Cp Vpa and its 
dual counterpart, 


rv _ _ _ kX 

^ p,v per ~ 


- p,v 


c, 


k,X per •> 


(5.1) 


as follows: 

The Weyl tensor vanishes identically for maximally symmetric backgrounds, 
Mink 4 , dS 4 or AdS 4 , which are solutions of the full non-linear Einstein equa¬ 
tions with cosmological constant A. Linearized gravity around the maximally 
symmetric solution, whose choice depends upon the sign of the A, exhibits a 
symmetry under the interchange of the electric and magnetic components of 
the corresponding Weyl tensor. At the level of equations, it amounts to ex¬ 
changing the role of the linearized field equations and the Bianchi identities. 
Two metric perturbations 8g liU and 8g lw are said to be dual to each other if 


Cpvpa ( g ) Cpvpcr ( g ) , 


(5.2) 


so that the electric components of the one Weyl tensor are the magnetic 
components of the other and vice-versa. Although these relations are highly 
non-local, they are naturally resolved by decomposing the gravitational per¬ 
turbations into two distinct classes of opposite parity, namely axial and polar 
perturbations. It turns out that 

5 ~polar = ^axial , ^axial = ^polar , ( 5 . 3 ) 


provided that the solutions of the corresponding effective Schrodinger prob¬ 
lems are interrelated as 

'hpol =- 'Lax , 'Lax = ‘Tr'bpol • (5.4) 

uj z 

Since the Schrodinger problems of the axial and polar perturbations are iden¬ 
tical, though the precise form of the effective potential depends on the cos¬ 
mological constant A, the relation among T and T enforces the two classes 


of perturbations to have the same frequency ui and the same boundary con¬ 
ditions at spatial infinity. 

Specializing the discussion to gravitational perturbations of AdS^ space- 
time, which is the subject of the present work, one finds that electric-magnetic 
duality has a holographic manifestation [2H - f2T] . Quite remarkably, the holo¬ 
graphic energy-momentum tensor T a b and the Cotton tensor C a b of the metric 
at the boundary of space-time follow from the electric and magnetic compo¬ 
nents of the bulk Weyl tensor (with respect to the radial ADM decomposition 
of the metric) as 


SirTab = - lim r 3 C arb r , 

(5.5) 

r—> oo 

C a b lilll T Carbr • 

(5.6) 


r—>-oo 


The Cotton tensor is traceless, symmetric and covariantly conserved, as for 
the energy-momentum tensor, and it measures the deviations of the boundary 
metric from conformal flatness due to the boundary conditions one imposes 
on the perturbations of the metric. The last equations are general and appli¬ 
cable to the full non-linear regime of AdS^ gravity, but when specialized to 
the linearized theory they give rise to the so called energy-momentum/Cotton 
tensor duality that is resolved in sectors as 

SrtT 1 " = Cf(5.7) 
ZkTT 11 = C". (5.8) 


5.2. Duality for Entanglement Entropy 

The extremal surface A extr , whose area provides the entanglement en¬ 
tropy, is defined with respect to a given bulk metric g for a given boundary 
region. Let us define A extr as the extremal surface taken with respect to the 
dual metric g, following (5.2), 


A extr (g) = A extr (g). 


(5.9) 


The corresponding entanglement entropy is taken to be 

J Area (i extr ) . (5.10) 

Since the unperturbed AdS± space is self dual, ^Ads = fl'AdS) the unper¬ 
turbed extremal surface is also self-dual. The variation of the dual entangle¬ 
ment entropy for axial and polar perturbations, which are dual to each other, 
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can be calculated in the same way as the variation of entanglement entropy. 
As a direct consequence of equations (5.3), we obtain 


<5S axial = (hS polar , (5.11) 

<5S' polar = <5S axial . (5.12) 


Similarly, one can define the variation of the dual modular Hamiltonian 
using the Cotton tensor of the original boundary metric, as 


5E 


^ I d ^ C 


<c 


(5.13) 


As direct consequence of equations (5.7) and (5.8), it turns out that 


<5£ ax ia i = h£ polar , (5.14) 

hE p oiar = <5£ axial . (5.15) 


It is true that 5S = 5E for all kind of perturbations, and, thus, it is also 
true that 

5S = 5E , (5.16) 

showing that the first law of thermodynamics for entanglement is self-dual 
under gravitational duality. 

Note in this context that the boundary metric is not identical to the dual 
boundary metric, since the duality transformation interchanges Dirichlet and 
Neumann boundary conditions for the metric (it only preserves the boundary 
conditions of the effective Schrodinger problems). Since the line element of 
the entangling curve remains inert to perturbations, the entangling curve is 
not identical to the dual entangling curve. Likewise, the polar cap region A 
is not the same as its dual counterpart A. These differences should be taken 
into account in the definitions of the dual thermodynamic quantities above. 


6. Conclusions and Discussion 

In this paper we studied the variation of the entanglement entropy and 
the expectation value of the modular Hamiltonian for small gravitational 
perturbations of AdS^ space-time satisfying general boundary conditions. 
We found that the first law of thermodynamics for entanglement is realized 
holographically through the Ryu-Takayanagi formula, provided that the line 
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element of the entangling curve on the boundary remains constant, hereby 
generalizing previous results on the subject. It should be viewed as addendum 
to the Ryu-Takayanagi prescription. In this context, we also noted that the 
perimeter of the entangling curve is the “area” appearing in the “area law” 
term for entanglement entropy. Thus, our demand that the line element of 
the entangling curve remains constant is equivalent to the statement that 
the “area law” term is inert to metric perturbations for all different kind of 
boundary conditions. In turn, this appears to be equivalent to using a certain 
renormalized form for the holographic entanglement entropy. 

Although our presentation is technically limited to metric perturbations 
and entangling curves that are axially symmetric, choosing polar cap regions 
with circular boundary, we expect it to generalize in the absence of axial 
symmetry. It will be interesting to see how exactly this happens. 

We also examined the implications of electric-magnetic duality of lin¬ 
earized gravity for the entanglement entropy and the associated first law of 
thermodynamics. We defined a dual entanglement entropy and a dual mod¬ 
ular Hamiltonian and found that they also satisfy the first law of thermody¬ 
namics for appropriately related dual regions A and A at the boundary of 
AdS^. These results too can be thought to provide an additional consistency 
check of the Ryu-Takayanagi conjecture. 

Typically, electric-magnetic duality connects the dual degrees of freedom 
of the theory in a non-trivial and non-local way. In the absence of excita¬ 
tions, it reduces to a trivial local relation, but this is not any more so in the 
presence of excitations. If one defines the entanglement entropy in a region 
by tracing out the degrees of freedom in the complement, the same degrees of 
freedom will correspond to a different region in the dual description. Thus, 
the dual entanglement entropy should be defined for a different entangling 
curve. Our results, showing that the holographic realization of the first law 
of thermodynamics is achieved only when the line element of the entangling 
curve remains constant, are in accord with the different boundary conditions 
satisfied by dual gravitational excitations so that the dual boundary metrics 
differ. 

It would be interesting to see if such considerations can also be applied 
to other backgrounds such as AdS^ black holes. Perturbations of black-holes 
exhibit a duality relation which is best described in terms of supersymmetric 
partner potentials for the effective Schrodinger problems of the axial and po¬ 
lar sectors. Although this duality has no direct interpretation in terms of the 
electric and magnetic components of the Weyl tensor in the bulk, as in AdS± 
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space-time, it has holographic manifestation as energy-momentum/Cotton 
tensor duality for the hydrodynamic modes of black holes [39]. We hope to 
address the related thermal effects of black holes elsewhere. 
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Appendix A. Useful relations in Poincare and Spherical Coordi¬ 
nates 

It is known that in Poincare coordinates, the extremal surface corre¬ 
sponding to a disk region with radius R centered at (x 0 ,y 0 ) is given by the 
“hemisphere”, 


(x - xaf + {y - y 0 f + z 2 = R 2 , 

r = T 0 . 


(A.l) 


Changing to global spherical coordinates, using the coordinate transforma¬ 
tion (2.13), the surface (A.l) takes the form, 


cos# = cos^m/l -|—- , 


t — to , 


where 9 0 and t 0 are given by 
cos 6* 0 = 


1 - R 2 + r 0 2 


cosf 0 = 


\j [(-R + r o) 2 + l] [{R ~ T o) 2 + l] 
_ 1 + R 2 - To 2 _ 

\J [(-R + r o ) 2 + l] [(-R - r o ) 2 + l] 


(A.2) 


(A.3) 

(A.4) 


32 









To express the variation of the modular Hamiltonian for the polar cap 
region in spherical coordinates, we needed the expression for the conformal 
Killing vector in the boundary, £, that leaves the entangling curve invariant. 
This is the limit of a conformal Killing vector in the bulk, £, which in Poincare 
coordinates and with the appropriate normalization takes the following form 

pnj, 

f ^ [(^ 2 - Z 2 - (t- r 0 ) 2 - (x - x 0 ) 2 - (y - y 0 ) 2 ) d T 

-2 (r - r 0 ) (. zd z + (x - x 0 ) d x J r (y- y 0 ) d y )\ . (A.5) 

In our analysis, the entangling curve is taken to be axially symmetric, as for 
the metric perturbations. This selects xq = yo = 0, so that the center of 
the disc is placed at the north pole. Furthermore, since the AdS^ metric is 
static, we may choose a specific instant of time to simplify the calculations. 
The choice r 0 = \/l + R 2 implies that 

tan 9 0 = R, (A. 6 ) 

cos to = 0 . (A. 7) 

At that instant, the conformal Killing vector £ in the bulk takes the form 


£ = - 
^ R 


R 2 — z 2 — [t — Vl + R 2 j — x 2 — y 2 j d T 

-2 (r - Vl + R 2 ^j (zd z + xd x + yd y ) 


(A. 8 ) 


Converting to spherical coordinates by the transformation (|2.13j) , we hnd 

2irVl + R 2 




R 


r cos 6 sin t 

Vl + r 2 


VTTr 2 


d t 


u , 2 a . Vl + r 2 sind cost 

— V 1 + r z cos t) cos to r H- Oq 


. (A.9) 


which can be rewritten as 




2tt 

sin 9n 


r cos 9 cos (t — t 0 ) 

Vl + f 2 


- cos 6*0 d t 


A ——2 a ■ u Vl + r 2 sin 9 sin (t — t 0 ) „ 

T v 1 + r 2 cos 9 sm (t — t 0 ) o r - o g 


(A. 10) 
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Then, the boundary limit of this expression provides the desired conformal 
Killing vector £, which reads 

2n 

C = ——— [(cos (t — t 0 ) cos 6 — cos 6* 0 ) d t — sin (f — t 0 ) sin^c^] (A. 11) 
sin c/q 

and is used to compute the variation of the modular Hamiltonian in section 
12.51 

Appendix B. The Modular Flow for a Disc in Minkowski Space- 
time 

It is not known how to express the modular Hamiltonian for a general 
region and state. This difficulty stems from the fact that, in general, the 
modular Hamiltonian is a non-local operator. Exceptionally, in some cases, 
the modular Hamiltonian generates a geometric flow and can be expressed 
as a local operator. Here, we outline how this is actually realized for a disk 
region in Minkowski space-time. 

The modular Hamiltonian defines a symmetry of the system. More specif¬ 
ically, the symmetry group is provided by the unitary operators U (s) = e~ lHs 
and it is called the modular group, 

Tr (pU (s) OU (s)) = Tr ( pp is Op ~ is ) = Tr ( p~ is pp is O ) = Tr (pO ). (B.l) 

We define O (s) = U (s) OU (— s). When the modular Hamiltonian is not a 
local operator, this flow is non-local. 

If one extends the modular group to imaginary parameters, there will be 
a periodicity relation in imaginary time, 

Tr (pOi U i 0 2 ) = Tr (pU (?) OiU f-r ! 0 2 ) = Tr (pp~ 1 0 1 p0 2 ) = Tr (p0 2 0i) . 

(B.2) 

Thus, the state p is thermal with respect to the time evolution of the modular 
symmetry with temperature T = 1. 

Let us consider Minkowski space-time in coordinates X ,L . The Rindler 
space 1Z is defined as the causal development of the half plane X 1 < 0. The 
causal development of a region A is the set of space-time points whose causal 
space-time curves necessarily cross region A. The modular transformations 
map the local operators O (A) in Rindler space 1Z to themselves. More 
specifically, the modular flow acts on the coordinates as follows, 

X ± (s) = X ± e ±2 ' ns , X i (s) = X\ (B.3) 
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where X ± = X 1 ± X° and the Latin indices i run from 2 to d — 1. 

The modular Hamiltonian for 1Z can be calculated explicitly. Transform¬ 
ing to the usual Rindler coordinates X ± = ze ±T / R , the metric becomes 

A 

ds 2 = ——dr 2 + dz 2 + dX i dX i , (B.4) 

R z 

corresponding to a thermal state with temperature T = 1/2itR. Thus, the 
density matrix can be expressed as 

^— c 1'kRH t 

PR — T r e -2irRHr (B-5) 

with modular Hamiltonian 


H n = 2 ttRH t + log Tie ~ 2nRHT . (B.6) 


The result for Rindler space can be subsequently used to obtain the mod¬ 
ular Hamiltonian for a disc D of radius R. The conformal transformation 


x p = 


X p - X V X V C P 


+ 


C p 


1 - 2X k C k + X p XPC a C a 2C\C X ’ 


(B.7) 


with 

C"= (°, T,°,,,,, 0 ) (B.8) 

yields 

r] lw dX ll dX 1 ' = 9 2 rj ll JX‘dx 1 ' (B.9) 

with 

Q=l- 2 X K C* + X p X p C a C a . (B.10) 

Letting r = yj(x 1 ) 2 + ... + ( x d_1 ) 2 and t = x°, we see that this conformal 
transformation maps the half plane A" 1 < 0 to the disk 1). r < R. and the 
Rindler space X ± < 0 to the causal development of the disk T>, x ± < 11. 
where x ± — r ± t. For later use, we express the new radial coordinate as 


r =U x ' xi+ { R -^§-)*- < b '°) 

Also, the inverse coordinate transformation is 

X p = O (x» + 2x u x u C p ) - , (B.12) 


35 



where 


n = ( - + c k x k + C\C x x v x v 


-i 


(B.13) 


We want to determine the form of the modular flow in the new coordinate 
system. For this purpose, we first make use of equations (B.3) and (B.12) to 
find the flow of the conformal factor fl, 

n(s) = 1 _Wd + -^ 


R 


4 R 2 


= Q — — [(cosh 2ns — 1) X 1 + sinh27rsA°] 


= ni- 

O 


1 R 2 + r 2 -t 2 


R 2 


+ - cosh 27 TS 


t 2 + R 2 - r 2 
R 2 


— sinh 2ns— 
R 


(B.14) 


= — (R cosh 7 ts — x + sinh ns) (R cosh ns + x sinh ns) 

RZ \ 


The radial coordinate r, which is given by equation (B.ll), also flows through 
its dependence on the conformal factor, 


r (s) = 


n(s) 


' X i X i + R- 


X u X l 
4 R 


Q 


n(s) 


(B.15) 


Similarly, the modular flow of the new time coordinate can be expressed as 
A" 0 (s) cosh 27 tsA° + sinh 27TSA 1 


t (s) = 


n(s) 

o 


R 


fi(s) 

fcoeh 2,rA-i sinh 2^ ^ 

R 2 


ft(s) V 
which in turn implies that 


Q (s) , _ ( t 1 

x (s) = r ± /t cosh27rs—-sinh '2ns 


(B.16) 


R 2 


t 2 + R 2 - r 2 


Q 




R 


R 2 


= — (f? cosh ns ± x T sinh ns) (x ± cosh ns R sinh ns) . 
R 


(B. 17) 


Finally, we obtain 

± x ± cosh 7 ts =F R sinh 7rs (R + x ± ) — (R~x ± )e :f27rs . . 

OC (s) — R— i i i ~ tt , j r ~ 

R cosh 7rs =F x ± smh ns (R + x±) + (R — aA) e =F2?rs 
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The modular flow in the new coordinates yields by differentiation 

(B.19) 
(B.20) 
(B.21) 

clarifying the action of the modular Hamiltonian on local operators in the 
new coordinates. Restricting attention to the time slice t — 0, one recovers 
the usual expression for the modular Hamiltonian for a disk of radius R jT8j , 
namely 

r f > 2 — r 2 

H d = 2tt / d^x —— T 00 . (B.22) 

J 2 R 


dx ± (s) 


ds 

dr (s) 
ds 

dt (s) 


= ±7T 


R 2 — x 


±2 


s=0 


R 


= -2vr 


s=0 


rt 

h ’ 


ds 


R 2 ~r 2 - t 2 


= 7T- 


s=0 


R 


Appendix C. Coordinate Dependence of Subleading Terms in <5A® xtr , 

S A extr 

The individual contributions to the variation of the area, <5A® xtr and 
(5A^ xtr , which appear in section |4j are coordinate dependent, unlike their 
sum which is coordinate independent when the line element of the entan¬ 
gling curve remains constant. Furthermore, the term <5A^ xtr only vanishes 
for Dirichlet boundary conditions provided that one works in coordinate sys¬ 
tems, unlike ours, where specific metric elements (to be made more precise 
later) contain no contributions that fall like 1/r close to the boundary. These 
statements are substantiated below by performing a small perturbative co¬ 
ordinate transformation altering the 1/r asymptotics of the metric elements, 


r = \Jr' 2 + e ( t , 6) r'. (C.l) 

The area of the extremal surface in AdS^ is provided by equation ( |2.21[ ). 
Writing the result in terms of the r' coordinate, we have 


^extr 



+ £ (to, 6 (r)) r' sin do - 1 j 


27r lim 

r'—> oo 



(C.2) 
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The change of coordinate introduces a change in the area of the minimal 
surface in AdS 4 , which functionally resembles <5A® xtr , 

SA' g extr = 7 re (t 0 , 9 0 ) sin 9 0 . (C.3) 

This seems to be odd at first sight, because the area should be independent 
of the choice of coordinates, but A extr is infinite and the result depends on 
way one is taking the limit. Ultimately we are interested in finding a result 
that is independent of the limiting process. 

On the other hand, the metric element g w of AdS 4 space in the r' coor¬ 
dinate is expressed as 

5W = { r ’ 2 + £ (*>#) r> ) sin261 • (C.4) 

Demanding that the length of the entangling curve remains invariant, we 
obtain 

SO (t 0 ) = -^—.£ (to A) tan 9 0 , (C.5) 

zr 

which, in turn, gives a non-vanishing contribution to the area of the minimal 
surface, which resembles 8 Ab CxtT , 

8 A' b extT = —ns (t 0 , 60 ) sin 9 0 . (C.6) 


Thus, demanding that the line element of the entangling curve remains 
constant, we arrive at the following result, 


8 A' 


+ $A' h 


= 0 . 


(C.7) 


This provides the desired coordinate independent prescription for taking the 
limit. The variation of the area of the minimal surface in AdS 4 should be 
zero under changes of coordinates. Likewise, for perturbations of AdS 4 space- 
time, the combined variation 8 A extl + 5A ^ tr , which is not zero, does not 
depend on the choice of coordinates, unlike the individual terms SA extT and 
8 Al xtr that depend upon it. 

As an illustrative example of the effect that coordinate transformations 
may have to the components of the variation of the area of the minimal 
surface in the presence of gravitational perturbations of AdS 4 space-time, 
we consider the case of polar perturbations. Using the change of coordinate 


(C.l) with the choice 


£ (t, 9) = e~ lut J 0 (-/(/ + 1) - oU ) P t (cos 9), 


(C.8) 
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we find 


5A ,extI = —7ie~ iujt ° lim 

^ r'—> oo 


r 'j\ + - Jo l (l + 1) ) sin 0 Q Pi (cos 6 0 ) 


■cotd 0 ^ Y~ J o ( P 1 +1 (c°s0o) - ^-1 (cos 9 0 )) 


(C.9) 


and 

5^4 / f) extr = ire~ lut ° lim r' J\ sin 9 0 Pi (cos 9 0 ). (C.10) 

r f —¥ oo 

One can readily verify that 8 A' g extI + <5^4^ extl = 8 A g extT + 8 A b extI , using the 
results of section [4j as advertised. 

Choosing Dirichlet boundary conditions for the metric, Ji = 0, we find 
that (5^4^, extl vanishes, whereas 8A' g extl is finite, equal to 4tt5E. Note at this 
end that the asymptotic behavior of g w in r 1 coordinate is 

9vp = r ' 2 sin 2 (l + O ^ . (C.ll) 


Summarizing, the term 8 Al xtT can only be neglected for Dirichlet bound¬ 
ary conditions and in coordinate systems where the metric elements cor¬ 
responding to elementary lengths parallel to the entangling surface do not 
contain subleading terms that fall like 1 jr close to the boundary. This is 
precisely the content of the calculations reported in earlier works mm, 
using Dirichlet boundary conditions and writing the metric in Fefferman- 
Graharn coordinates which obey the restrictions stated above. If any one of 
these restrictions do not hold, the variance 8 A™ U will become relevant in the 
calculation and the need for invariance of the line element of the entangling 
surface will come into play, as in our calculations. As emphasized above, the 
individual terms 8 A g extr and 8 A b extT contain gauge artifacts, but not their 
sum, which is coordinate independent. 
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